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Abstract: In this communication we introduce a new model which represents the 
interaction between an atom and two fields injected simultaneously within a cavity in- 
^ eluding the nonlinear couplers. By using the canonical transformation the model can 

^ ' be regarded as a generalization of several well known models. We calculate and discuss 

O^. entanglement between the tripartite system of one atom and the two cavity modes. For a 

^ short interaction time, similarities between the behavior based on our solution compared 

' with the other simulation based on a numerical linear algebra solution of the original 

m : 

CN . Hamiltonian with truncated Fock bases for each mode, is shown. For a specific value of 

O ' the Kerr-like medium defined in this paper, we find that the entanglement, as measured 

^ I by concurrence, may terminate abruptly in a finite time. 

OO 
O 

PACS: 32.80.-t; 42.50. Ct; 03.65.Ud; 03.65.Yz. 

- - . 1 Introduction 

Recently, the experimental generation of atom-field entangled states in the mesoscopic 
field regime has been reported for fields with an average photon number of a few tens 
pp. Several entanglement measures have been studied, for example the von Neumann 
reduced entropy |2], the relative entropy of entanglement [3]. Several authors proposed 
physically motivated postulates to characterize entanglement measures [2]-[5]. Although 
these postulates vary from one author to another in the details, however they have in 
common that they are based on the concepts of the operational formulation of quantum 
mechanics [6]. The entanglement properties have been reported in many different cases 
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On the other hand, various generahzations of the atom-field interaction models have 
been considered. One such generalization is to inject two fields simultaneously within 
a high-Q two-mode bichromatic cavity. In this case the atom interacts with each field 
individually as well as both fields [0]. In the direction of nonlinear generalization of the 
atom-field interaction the infiuence of nonlinear coupling of a cavity mode with a Kerr- 
like medium on the decay of the excited state of the atom was reported [TOl Hi]. The 
nonlinearity can make the dynamics more intricate, for example with respect to switching, 
modulation and frequency selection of radiation in optical communication networks [121 
[131 [Hj. Ill addition, the presence of a nonlinear medium is of particular interest as atom- 
field dynamics is significantly affected. 

This encouraged and stimulated us to investigate a new version of the atom-field 
interaction model by including the nonlinear couplers. The main goal of this work is 
to investigate the effect of the nonlinear medium on the phenomenon of entanglement 
when two fields simultaneously injected within a high-Q cavity. We introduce a way to 
study tripartite entanglement which is an important issue, by combining two fields and 
an atom in a double cavity geometry with a Kerr nonlinearity. We also discuss how the 
nonlinear effects increase in importance compared with increasing atom-field coupling. We 
show that nonlinear considerations can play a major role in determining the maximum 
entangled states of the system and generating C-Note gate. Particularly, we employ linear 
entropy and concurrence to quantify coherence and entanglement, respectively. Beyond 
fundamental investigations, these results may be of benefit in the characterization of atom- 
field parameters, and controlling these features may allow for new control techniques for 
single and multiple qubit coupling. 

The paper is organized as follows: In section 2 we describe the model and the analyti- 
cal approach which based on the canonical transformation. In Sec. 3 we briefiy comment 
on the collapse-revival phenomena. In Sec. 4 we identify the different regimes of en- 
tanglement due to the linear entropy and concurrence. Sec. 5 presents some concluding 
remarks. 

2 The model and its solution 

We devote this section to introduce our model taking into consideration the effect of a 
Kerr-like medium (nonlinear couplers). For this reason let us make our starting point the 
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Hamiltonian 

H = Hf + H, + H,i, (1) 

where Hf comprises a field-field interaction (parametric down conversion model) in the 
presence of a nonlinear Kerr medium given by 

2 

H{ = ^ H ^(^Uidldi + Xja^r^i) + XCi\aid\a2 + h\ {a\d2 + a\a^ , (2) 

i=l 

where Xi?^ = 1)2 and x are related to the cubic susceptibility of the medium, such that 
Xi,i = 1,2 represent the self-action for each mode, while x is related to the cross-action 
processes, respectively [I6] . The other two parts of the Hamiltonian consists of the atom- 
field interaction and the free atomic system that is 

2 

i^a = -^^^^^ -^'^f ^ ^-^^ (aicr+ + a|a_ j . (3) 

i=l 

The Hamiltonian ([2]) can be regarded as a generalization of models considered earlier in 
the absence of the Kerr-like medium, see for example [T7j, which concentrate on studying 
entropy squeezing as well the degree of entanglement. In order to consider the statistical 
properties of the present system we have to find the dynamical operators either by solving 
the equations of motion in the Heisenberg picture, or by finding the wave function in the 
Schrodinger picture. 

The latter case will be adopted. To reach our goal let us first introduce the canonical 
transformation 

cos 6 sin 6 
- sin 9 cos 9 




) 






{ b2 j 



(4) 



with the properties 



a[di + 020-2 = ^1^1 + ^2^2, (5) 



where [04,0.]] = [bi,b^j] = ^ij- Since the total energy for the system before the rotation 
is equal to that after rotation and hence the transformation is invariant. Thus, with a 
particular choice of the angle 9, the above transformation would enable us to remove the 
evanescent waves term from the Hamiltonian ([2]). This can be achieved if we take 

^ = -tan-if^^V (6) 

In this case the Hamiltonian ([2]) reduces to 

2 2 
H = Y.n [{n}>% + xbfbi) + /i. {b^a+ + b\a^)] + % b%b]h + y^^^, (7) 

1=1 i¥=j=i 
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where we have assumed the cross- action couphng parameter is equal to twice the self- 
action parameter of each mode, so that Xi — X2 — X — (codirectional coupler case). 
Further, we have also defined 



Qi = {ui cos^ 9 + u;2 sin^ 61 - A sin 26) % Q2 = {<^2 cos^ 6 + u;i sin^ 61 A sin 26) ^ , 
and 

fjL^ — Ai cos 6 — X2 sin 6, — A2 cos ^ -I- Ai sin 6. 

Finally, let us adjust the coupling parameter A to take the form 

'(jJ2 - OJl 



(8) 
(9) 



A — A1A2 . 2 ^2 ) ' 

V A2 - Ai / 

from which the interaction Hamiltonian V^{t) can be written thus 

V\t) = Ji {h2&+ ® exp (^-i[A + 2x{h\W + 1% - l]t^ + c.c) , 



(10) 



(11) 



where /i — -\/A^ + A2 and A = — o^o- Having obtained the interaction picture, we are 
therefore in a position to find the explicit solution of the wave function and consequently 
the density matrix. After some calculations the density matrix of the system can be 
written as 



p{t) = u{t)pm\t), 



(12) 



where U{t) is a 2 x 2 matrix representing the time evolution operator; its elements are 



COS d^^^rn^t + ^— + sm 5^^^^^t 

sin^Jl, 



exp( i^mi,m2^)i 



U22{t) 



COS d^^^^^t - I— sm 5^^^„J. 

^m\,m2 



cxp(ir' _i). 



(13) 



'mi,m2 



m\,m2 



^-^ nil, 1712 



f{m2 + 1) 



' m\ ,m2 



mi,m2 



p?m2 



= A + 2x(mi +m2 - 1), F; 



mi,m2 



A + 2x(mi + m2-2). 



(14) 
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3 Atomic inversion 



For applications in real systems, one can see the atomic inversion of the two-level system 
is of particular interest. Therefore, in the present communication we discuss the time 
dependence of the two-level-system observable {az{t)), for different values of the Kerr-like 
medium. The usual coherent state is used as initial conditions for the fields. As might 
be expected, the behavior of the two-level system changes dramatically depending on the 
value of the non-linear medium. 

1 
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Figure 1: Atomic inversion as a function of the scaled time Ait. The parameters are 
A2/A1 = 0.01,cJi/Ai = 0.2,cj2/Ai = 0.1, A/Ai = 0,x/Ai = 0, and = = 10. The 
initial state is \ai,a2]e), where = |aj|e*'^, x/^i = 0.001 (sohd curve) and x/Ai = 0.1 
for (dots). 

In Fig. (pP) we have plotted the atomic inversion as a function of the scaled time Ait 
for different values of the Kerr-like medium and the phase of coherent states is taken to 
be zero i.e. = 0. It has been observed that for a small value of the Kerr-like medium 
such as x/^i = 0.001, the system shows period of collapse after onset of the interaction. 
This is followed with a long period of the revival. Here we may point out that comparing 
with the first period of collapse the amplitude of the oscillations for the second period 
is decreased, see fig.(l) (solid curve). This means that as the time of the interaction 
increases the amplitude of the oscillations decreases. Increasing the value of the Kerr-like 
medium leads to increase in the period of collapses with decreasing in the size of the 
amplitude of each period, see fig.(l) (dot curve). It is also noted that as a special case of 
the model when A2 = X2 = O5 general behavior of the atomic inversion coincides with 
that of the well known JCM in presence of the Kerr-like medium [18|. 

We have verified the behavior obtained in the numerical simulations based on the final 
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equations and the numerical solutions of the original Hamiltonian i.e. one simulation 
simply based on the final equations and the other simulation based on a numerical linear 
algebra solution of the original Hamiltonian with truncated Fock bases for each mode 
(see Fig. [2]). Provided the ratio A2/A1 is made small enough, we have excellent agreement 
for the short time scales. In the meantime, the value A2/A1 = 0.01 gives quantitative 
agreement over some periods of oscillation. 




-0 .75 

-1 



Figure 2: Atomic inversion as a function of the scaled time Ait. The parameters are 
A2/A1 = 0.01,tui/Ai = 0.2, x/Ai = 0.01 U2/\i = 0.1, A/Ai = 0, x/Ai = 0.01. The initial 
state is |5,6;e), where dotted curve obtained using the final equations and solid curve 
obtained using linear algebra solution of the original Hamiltonian. 

The revival time for the present system is calculated to give us 

t (15) 
A {^n — ^n-l) 

where fin is the Rabi frequency, such that fin = m2 integer, 

2n7r f 1 



+ -\/(2x'(2n-2) + A')2+7?%j, (16) 

where x' = x/'^ii 'P' = 'P'/^i, ^' = A/Ai- In the following section we turn our attention 
to discuss the effect of the Kerr-like medium on the phenomenon of entanglement. 



4 Entanglement 

In this section we shall concentrate on the discussion of the entanglement where the total 
state vector can not be written precisely as the product of a time-dependent atomic and 
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field component vector. 



4.1 Coherence loss 

Distillable entanglement is a critical resource for quantum information. In our case we are 
considering transduction of quantum information between fields and atoms so the degree 
of entanglement between the atom and the field is important to assess its use in quantum 
information. Furthermore there are three entities: two field modes and one atom, and the 
entanglement between any single entity and the other two can be used so we assess the 
inherent resource in our system by calculating the tripartite entanglement in the combined 
system. 

Here we use the idempotency defect [19] , defined by linear entropy, as a measure of the 
degree of mixture for a state Pa{t)- This version of entropy makes the problem tractable 
[20] and therefore we use this idempotency defect as a measure of coherence loss. This is 
given by 

4"^ = Tr[p,(t)(l-p,(t))], (17) 
where 8j:"'^ has a zero value for a pure state and 1 for a completely mixed state. 




Alt Alt 



Figure 3: Linear entropy as a function of the scaled time Ait. The parameters are A2/A1 = 
0.1,ci;i/Ai = 0.2,co'2/Ai = 0.1, A/Ai = 0,x/Ai = 0, and ni = n2 = 10. The initial state is 
7f(|e) + Ig)) ® |tti,a2), where (a) x/Ai = 0.001 and (b) x/Ai = 0.01. 

In Fig. ([3]) we have plotted the idempotency defect as a function of the scaled time 
for different values of the Kerr-like medium. It is shown that the asymptotic value of the 
linear entropy is obtained when the time is increased. Also the value of Sj;""^ increases 
as the time increases in addition to the appearance of irregular fluctuations at discrete 
period of time for a large value of the Kerr-like medium, see fig. (3b). Of course, there are 
some differences between the two cases (small and large values of the Kerr-like medium) 
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in the amplitudes but the general behavior is the same. This may be thought to arise 
from the asymptotic limits which have been observed in both Figs. [3]^a,b). 

4.2 Tripartite Quantum States 

Now we would like to discuss an example of entanglement between the two cavity fields 
and atom. In the case of pure state, if the density matrix obtained from the partial 
trace over other subsystems is not pure the state is entangled. Consequently, for the pure 
state of a bipartite system, entropy of the density matrix associated with either of the two 
subsystems is a good measure of entanglement. In the mixed state case, the entanglement 
can be quantified by the quantity C{p) which is known in literature as concurrence. Quite 
recently, some approaches have been reported for the determination of entanglement in 
experiment [211 [221 [23|, [2l] . The most remarkable are the new formulation of concurrence 
in terms of copies of the state which led to the first direct experimental evaluation of 
entanglement and some analogous contributions [2S] to multipartite concurrence. 

For the density matrix p(t), which represents the state of a bipartite system, concur- 
rence is defined as [26] 

C(p) =max{0,Ai- Aa- As- A4}, (18) 

where the A, are the non- negative eigenvalues, in decreasing order (Ai > A2 > A3 > A4), 
of the Hermitian matrix 




Here, p* represents the complex conjugate of the density matrix p when it is expressed 
in a fixed basis and ay represents the Pauli matrix in the same basis. The function C (p) 
ranges from for a separable state to 1 for a maximum entanglement. 

We first investigate the quantum correlation between the atom and cavity modes. If we 
consider the initial state of the system is given by p(0) = 7|0, 1; e) (0, 1; e|+/3|0, 1; (7)(0, 1; 
7 + /5 = 1, which means that the atom starts from a mixed state and the field starts from 
|0, 1) state i.e. the vacuum for the first mode and one photon in the other mode. If we 
deal with the two cavity modes as system A, and the atom as system B, then p{t) in 
equation fll2p can be thought of as the density operator of a two-qubit mixed state. In 
the basis |1) = |01) O |e), |2) = |01) ® \g), |3) = |00) O |e), |4) = |02) \g), the density 
matrix can be written as [27j 
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Pit) 



(20) 



/ \ 

P22 P23 

P32 P33 

\ P41 P44 / 

where p^^- = {i\p{t)\j). 

The exphcit expression of the concurrence describing the entanglement between the 
system A and system B can be found using equations (fTSjl and (120|) as 



C(p) = 2max{0, IP23I - VP11P44, IPmI - VP22P33} 



(21) 




Figure 4: Concurrence as a function of the scaled time Ait. The parameters are A2/A1 = 
0.1,Wi/Ai = 0.2,cj2/Ai = 0.1, A/Ai = 0,x/Ai = 0. The initial state is p(0) = |(|e)(e| + 
|^)(^|) (g) |0, 1)(0, 1|, where (a) x/Ai = 0.001 and (b) x/Ai = 0.01. 

The graphs of the concurrence as a function of time for the two dynamical regimes 
are displayed in Fig. 4. Keeping the value of the Kerr-like medium small enough, the 
asymptotic value of the concurrence is not null, since the global system evolves to a 
classically correlated state. But once the x/^i is increased we see that the concurrence 
vanishes in the asymptotic limit. With increasing x/^i further, we note that the decay of 
the concurrence is more rapid than the corresponding decay in the weak nonlinear regime. 

If the initial state has been considered such as 7 = 1, then the sudden death of 
entanglement will be obtained if P11P44 = 0. This means that the sudden death time is 
given by 

1 ^ f 2 (A' - 2x'f \ , , 

Aitrf = cos"^ ^ ^ . 22 

2x'-A' l^(A'-2x')'-8(l + A'2^); ^ ^ 

In the limit that the fields decouple from the atom, it is shown that, one may just entan- 
gling the cavity fields along the lines which has been considered in Ref. j25] . 
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In quantum computation operations are performed by means of single-qubit and 
multiple-qubit quantum logic gates. In what follows we consider universal quantum logic 
gate based on two electromagnetic field modes of a cavity. From equation (fT2l) . we can 
get various time evolutions of the present system. If the initial state is taken to be |e), 
and we set 

X = yv/4^, (23) 

then, to obtain the C-NOT gate, we need the unitary operation with x = ^\/4 — /J^, 
A = 0, and interaction time Ait = 2mT, {n = 1, 2, 3, ...) which gives 



|0,0)/K 


^ |e)a|0,0)/h- 


^ |e)a|0 


J0,1)/^ 


-I^7)a|0,l)/H 


- I^7)a|0 



|e)a|0,0)/ ^ 

\g)a\0,l)f ^ \g),\0,l)f^\g),\0,l)f^\g)a\OA)f, (24) 

i.e. using a specific values of the Kerr-like medium in the present model, one can be able 
to implement a C-NOT gate. 



5 Conclusion 

In this paper we have discussed a new Hamiltonian class which describes the the inter- 
action between a single two-level atom and bimodal cavity field taking into account an 
optical Kerr nonlinearity. We present an analytically solution and numerical investigation 
of the atomic inversion, coherence loss and entanglement. Entanglement of the tripartite 
system of one atom and the two cavity modes has been discussed and the sudden death 
of entanglement is shown. Also, we proposed a scheme for quantum computing (C-NOT 
gate), which is realized by a nonlinear interaction in the QED cavity. For our model to be 
useful for experiments, we would need to include cavity losses, fluorescence, and atomic 
motion. These extensions to our model could be accomplished by standard master equa- 
tion methods, and such an analysis is reserved for future work. Our purpose here has 
been to develop a Hamiltonian that includes many of the iconic Hamiltonians used in 
quantum optics as special cases so that three coupled systems are used, which opens up 
investigations from bipartite to tripartite entanglement. Cavity quantum electrodynamics 
with one atom and two fields, including a nonlinear interaction, creates rich, interesting 
dynamics as shown here. 
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